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The mode-mode entanglement between trapped ions and cavity fields is investigated in the disper- 
sive regime. We show how a simple initial preparation of Gaussian coherent states and a postselection 
may be used to generate motional non-local mesoscopic states (NLMS) involving ions in different 
traps. We also present a study of the entanglement induced by dynamical Stark-shifts considering a 
cluster of N-trapped ions. In this case, all entanglement is due to the dependence of the Stark-shifts 
on the ions' state of motion manifested as a cross-Kerr interaction between each ion and the field. 
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When thinking of practical implementation of quantum computers and protocols, two types of information carriers 
have been widely considered. The first type is composed of two-level systems (qubits) such as quantum dots [l[ or 
superconducting devices Q , and the second type is composed of continuous variable CV systems like the continuous 
quadrature amplitudes of electromagnetic fields Q. These two types of quantum information carriers are usually 
thought as unit blocks (flying and stationary qubits) for distributed quantum computation Q. For long distance 
communication, photons are more appropriate to carry information between stationary qubits and are labeled flying 
qubits. It is then clear that in order to fully explore the distributed quantum computation content, it is important to 
' investigate the interaction between CV systems and qubits. In this context, trapped ions interacting with quantized 
£vq , electromagnetic fields d, H, 0, § constitutes a promising physical system for it comprises the interaction of a finite 
system (internal electronic levels) with two bosonic modes namely the cavity field and the vibrational center-of-mass 
motion. Distributed quantum computation relies heavily on entanglement which is one of the key ingredients for 
the achievement of quantum information tasks. For instance, when supplied with a pair of maximally entangled 
states, two parties may achieve perfect teleportation @ and dense coding [l(| using local operations and classical 
communications (LOCC). Such applications require entanglement as a resource and clearly show the importance of 
' studying it in different physical systems for actual implementation. 
Oh! In this work, we study the entanglement properties of a cluster of two-level trapped ions interacting with a quantized 
h-^ ' electromagnetic field. This multipartite system has been first considered by Buzek et al in Q where they explore a 
special resonance condition allowing the system to be effectively described as a collective Tavis-Cummings model [ll| . 
Here, we consider the opposite situation where the cavity and the internal electronic levels are kept far off resonance 
q-( with each other so that the field interacts only dispersively with the ions. This paper is organized as follows. In 
• • Section |TT] we describe the dispersive interaction of N-trapped ions with a quantized electromagnetic field. In Section 
IIIII we show that depending on the initial state of cluster and field, the cavity induced atomic dipole-dipole couplings 
may lead to the generation of non-local mesoscopic states (NLMS) [HI, [l3l. Jbj - These states are known to violate 
Bell inequalities for suitable choices of sampling points in phase space [15j . and in our generation protocol, they 
are naturally protected from destructive effects of cavity losses. The NLMS generation is achieved after a simple 
measurement of the internal state of the ions (postselection). Section ITVl is dedicated to the study of mode-mode 
entanglement induced by motion-dependent Stark-shifts (l6| arising from the dispersive interaction with the field. 
In this situation, we analyze the dynamics of entanglement between the elements in the cluster of ions or between 
field and cluster, and compare these bipartite quantities to the multipartite global entanglement measure [l7| whose 
evaluation involves all ions and the cavity field. The ideas and effects discussed here in terms of trapped ions could, 
in principle, be extended to neutral atoms trapped in optical lattices [la ]. The standing wave laser beams possess 
the periodic structure able to trap individual atoms the way it is needed to the applications shown in this paper. 
Of course, an additional experimental challenge would be to put such atoms in interaction with a quantized field, 
typically a cavity mode. Section [V] is devoted to some conclusions and final comments. 



2 



II. MODEL 



We consider a cluster of N-trapped ions placed in a resonator which sustains a particular mode of electromagnetic 
field. The whole system is treated quantum mechanically and its Hamiltonian is given by(H =1)0] 

N N „ z N 

H = ^ Via\ai + u c Pb + Y UJi ~2 + ^ giC0S [ 7?i W + &i)]{°'tt > + ^i^)> (!) 

2—1 2—1 i — 1 

where denotes the annihilation (creation) operator for the vibrational mode of center-of-mass motion of ion i 

(frequency i/j), bi(b\) denotes the annihilation (creation) operator for cavity field mode (frequency oj c ), &i's are the 
usual Pauli atomic operators for the two-level trapped ion i, is the field- ion i coupling constant, and rji = 2nAxi/X 
is the Lamb-Dicke parameter, Axi being the rms fluctuation of the position of ion i in its lowest trap eigenstate, and 
A the wavelength of the cavity field. Also, we are considering that the traps are sufficiently far apart so that one can 
neglect the Coulomb repulsion between the ions. 

In what follows, we will be interested in the cross-Kerr coupling of the bosonic field and motion operators. This is 
derived in detail in [l6( for the case N = 1 where the dispersive interaction of ion-field has been considered. We are 
now going to derive the correspondent effective Hamiltonian accounting for the dispersive interaction of the cavity 
field and the cluster. This may be achieved by a careful analysis of the quantum equations of motion and application 
of the high detuning condition A 4 = u>i — ui c ^> gi. We first note that the general Hamiltonian (fTJ is highly nonlinear 
because the function cos?y(aj + fij) contains powers of combinations of bosonic operators. It is shown in [l6l. [T^| that 

in the regime Aj -C v%, it is possible to keep just the terms proportional to powers of the number operator a\a,i so 
that the Hamiltonian (TTJ may be approximated by 

N N ^ z N 

H = Y, + u <$b + w *lf + 9if(atai)(o-tb + a^tf), (2) 

i—1 i—1 i—1 

where 

/(a+a) = e-" 2/2 : J (2r]V¥E) :, (3) 
and : Jo : is the normally ordered zeroth order Bessel function of the first kind [2(| ■ 

n 

N 

,,.h = A,/rtj> -4- n, ffnXnA^.W -I- ' 

dt 1 



The Heisenberg equation of motion for o~\b with 1< j < N when using the system Hamiltonian |T]) is given by 

ii & f b = A ^ b + Qjfiaja^afbP + J2 gifi&l^W^- ( 4 ) 



i=l 



It is now convenient to move to an appropriate reference frame by making — > a't e tuJjt , b 
dje~ lUjt . In this new frame, we may rewrite ^ as 

d N 

= 9jn^aj)afbbh- iA ^ + ^^/(aJaOara+e-^ 4 . (5) 

1=1 

In the limit A& 3> gk with 1 < k < N, we can integrate ([5]) to obtain in the original frame 

N 

J l—l 

Now, by substituting (J6J) in the Hamiltonian ((2]), we obtain the effective Hamiltonian pll ] 

H c ff = H + i?stark + H dip , (7) 

with the free part given by 

N N a z 

H = Y^l^i +u c Pb + Y l u} i-£-' ( 8 ) 

i=l i=l 
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the Stark-shifts given by 

JV 2 

ffsta rk = £ frA&frMCi + ^R^" - S^r^l (9) 

and the energy exchange term takes the form of a atomic dipole-dipole coupling given by 

JV 

H dip = J2 ^fi&laMiafaXaf&J + &-&f). (10) 

The effective Hamiltonian ([7} is valid in the regime Vk ^S> Afe ^> gk, with 1 < k < N. This effective Hamiltonian 
describes several interesting physical effects. We note in Q that the dynamical Stark-shifts depend upon the ions' 
center-of-mass motion via f{a\a,i). Despite the obvious interest this might attract from the point of view of funda- 
mental atomic and optical physics, it turns out to be an unusual way for entangling bosonic modes in a controlled way 
as we show later in this paper . We may also see that the dipole-dipole Hamiltonian also depends on the intensity of 
the vibrational motion of the ions via f(a\a,i). In the next section, we explore such Hamiltonian for generating NLMS 
involving the motional state of ions in different traps. 

It is worth noticing that the Hamiltonian ([7]) splits into ([9]) or (fTX)|) depending on the initial preparation of the 
internal state of the ions. This is a very useful feature when we think about the use of such a system for the 
implementation of quantum information tasks. For example, consider that all N ions are initially prepared in their 
ground state \g)® N ■ In this case, the effective Hamiltonian in interaction picture will be reduced to 

JV o 

Hlf N =-Y,y 2 (4am. (ID 

2=1 

This Hamiltonian describes a cross-Kerr interaction between the vibrational bosonic operators of any element of 
the cluster of ions with the bosonic field operators. The single ion case has been previously studied in [l6l ] where a 
quantum phase gate robust against spontaneous emission has been proposed. Again, it is noteworthy that the physical 
mechanism responsible for the appearance of the cross-Kerr interaction between motion and field is the dynamical 
atomic Stark-shifts induced by the quantized cavity field. As already said before, we will analyze the entanglement 
generated due to this coupling in Section HVl 



III. NLMS GENERATION PROTOCOL 



We now propose a scheme for generating entangled even and odd states of the form 

\* + ) = ^-(\C + ,C + ) + \C-,C-)) (12) 

or 

= j^(\C+,C-) + \C-,C+)), (13) 

where \C±) = (\a) ± | — a))/N±, with | ± a) being single-mode coherent states and N±^ + ^ + the normalization 
constants. The entangled states (fT2|) and (fT3|) are also commonly called nonlocal mesoscopic states NLMS in the 
case they involve two well spatially separate modes, and they are very suitable for the study of nonlocal properties 
of quantum mechanics. Actually, the states |C+) and |C_) are eigenstates of the parity operator e 47m+a , and this 
fact implies that the states (TT2")) and (fTB"]) involve the entanglement of different parities for the bosonic modes. This 
is the key ingredient for testing non-locality in these two- mode entangled states as shown in [l5j]. They propose the 
use of a adapted version of a simple Bell inequality that may be violated by the NLMS for carefully chosen points 
in phase space. Although, there are some other proposals in the literature for generating NLMS in the context of 
cavity quantum electrodynamics [Til . [23[ or trapped ions [HI , we consider here the system comprising trapped ions 
and quantized fields altogether, and this brings out a quite different way for accomplishing the generation of NLMS. 
This new scheme makes use of the atomic dipole-dipole coupling that depends on the intensity of the motion of the 
trapped ions (fTU|) . Due to the combined effect of the Rabi oscillations with the motional dependent coupling, the 
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whole system naturally evolves to a multipartite entangled state that, after postselection, collapses to the desired 
NLMS involving ions in different traps, as we are going to show next. Let us consider now the case iV = 2 in (0 and 
also that the electronic initial state of one ion has been prepared in excited state \e{) and the other in the ground state 
\g-z). Just like before, the effective Hamiltonian ([7]) reduces to a specific form depending on the internal electronic 
state preparation, and for the case now considered it reduces to 



ir|ff >ls3> = via\a\ + v 2 a\a 2 + oj c Pb + lj^ + w 2 y + (1 + bb r ) 
, flifl2 ( Ai + A 2 



2 2 
9l /2/"t» \--\ i 92 r2/~]~ \-+-- 

— f (aia^aZ <r x + {a 2 a 2 )a 2 r a 2 



(r££f) f(^)M^t^ + (i4) 

If we consider now the two ions to be identical, i.e, Ai = A 2 = A, the Hamiltonian (Tj"4"|) will be written in the 
interaction picture as 

£|e ffl >| 92 > = 9j^ f{ ^ i)f{ ^ 2){a + & - (15) 

We note that the Hamiltonian (| 1 5jl is independent of the state of the cavity field. It follows that applications relying 
upon (I15p will be robust against T = cavity losses and would work well even for ideal cavities sustaining a thermal 
field [|f7 

The NLMS generation protocol starts with the preparation of the system's initial state 

|V>(0)) = \e 1 ,a 1 ,g 2 ,a 2 ), (16) 

i.e, the vibrational degree of freedom of the ion 1(2) is prepared in the coherent state |ai( 2 )). We would like to point 
out that such initial preparation is quite feasible in the actual experimental status of laser manipulated trapped ions 
(2^ | . As a matter of fact, the experimental realization of vibrational coherent states for a 9 Be + ion has already been 
reported [13]. 

According to (fl"5]) the initial state (fT6|) evolves to 

= \e 1 ,g 2 )cos(Qt)\a 1 ,a 2 ) -i\g 1 ,e 2 )sin(nt)\a 1 ,a 2 ), (17) 

where 

n = ^f(&Ui)f(4&2). (is) 

We now assume the Lamb-Dicke regime (771,2 <C 1) where the motion of the ions is restrict to a spatial extent much 
smaller than the cavity field wavelength. For a further simplification, let us also take 771 = rj 2 — 77, and then the 
operator (fl"8|) may be approximated by 



ft ~ 9192 



(l- 77 2 - 7 7 2 [a t 1 a 1 -4a 2 ]) . (19) 



A 

In this approximation, the evolved state of the system (|17p will assume the form 

|V(t)> = \e 1 ,g 2 ){cos(uj ri t)[\<t>t, <t>+) + \<fc, <t> 2 )] - 1 sm(w^)[|^, <j>+) + \<f>+, <fq )]} 

+| 9 i, ea){cos(w,t)[|^, 4>t) + ^^2)} - * sin(w 77 t)[|0+ </>+) + \<fc, 4> 2 )]}, (20) 

where Lu n — <7i<? 2 (l — rj 2 )/ A, and 

I i V 2 9192 f , 1 - '? 2 8182 f , 

\<j) + ) '' ' " ' " " A ' 



j ' 2 



■») 8132 t , 1 j '? 81 82 



(21) 



' 2t 

with j = 1, 2. Now, taking ai — a 2 = a, and an interaction time v 9 \ S2 ti = 7r/2, the state (f20f takes the form 

|^(i/)) oc |*+){cos^|ei,5 2 ) ~ 7sin^|gi,e 2 )} + \$ + ){cosQ v \gi, e 2 ) -«sia^|ei,5 2 )}, (22) 



5 



(1 i 

with Or, = u ~V 7r - Now we will explain the post-selection procedure that will allow us to obtain the NLMS in the 



form (|12p and (jT3|) for the motion of the trapped ions. If the Lamb-Dicke parameter is carefully set to fulfill the 
condition 9, q — kir/2, i.e. r\ = + 2k), with the k integer or half-integer, the state ([2"2"| reduces to either 

mi)) cx |* + )| ei , 92 ) + |$ + )| ff x,e 2 ), (23) 

for the k integer, and 

\i/>(ti)) cx \* + )\g 1 ,e 2 ) + \<S> + )\e 1 ,g 2 ), (24) 

for the k half-integer. Then, after choosing k, a simple measurement on the internal state of the ions would lead 
to the generation of the states (fT2|) and (fT3|) . The experimental discrimination between the two internal electronic 
levels may be done using the very efficient electron shelving method [2§, 28]. Following our proposal, one could then 
generate NLMS involving the state of well separate trap ped ions and then perform the experimental investigation of 
its nonlocal properties using the proposals described in [131 . Finally, we would like to mention that there are other 
different mechanisms and regimes used to generate entanglement in this system [|| [2^| . 



IV. ENTANGLEMENT STUDY 



We now present a study of the entanglement properties of our model. When the cluster of ions is initially prepared 
in the ground electronic state \g)® N , the resulting Hamiltonian is the cross-Kerr interaction (|lip. It involves the 
coupling between the vibrational motion of any ion in the cluster and the cavity field. Just like before, we will be 
interested in the common Lamb-Dicke regime characterized by the condition r/i ^ 1. In this regime, the Hamiltonian 
(fTTj) will read 

N 

#lff *" = - J2 (25) 

i=l 

with Xi = 2n 2 gf /Aj. We remark that the natural frequencies of all oscillators have been shifted and included in the 
free Hamiltonian not shown [1 61 ] . 

In order to study the dynamics of entanglement in our model (|25p , we employ the linear entropy [30l | as a bipartite 
entanglement measure to analyze different partitions of the system, and the global entanglement (GE) [l~7| to discuss 
the multipartite entanglement between all elements of the cluster and the cavity field. The linear entropy and the 
global entanglement have been used quite extensively in the investigation of entanglement behavior in quantum phase 
transitions [3ll.l32l|. 

In what follows we investigate the entanglement dynamics for three types of initial states of physical interest. First 
we consider both the ions and the field in coherent states (case A) , then the ions in coherent states and the field in a 
squeezed vacuum state (case B), and finally the ions in squeezed vacuum states and the field in a coherent state (case 

A) Initially coherent states: let us consider the following initial state |^(0)) = <8>™ = i|oy) eg) |ao). Since the interaction 
is diagonal, the time-evolved state \^(t)) = U(t)\^(0)) can be easily calculated and the reduced density operators for 
one of the ions (fc-th, k = 1, ..N) and the field (tracing over all the other degrees of freedom) can be exactly calculated. 
The linear entropy for the a-th subsystem can be calculated by means of the expression Sl(Po) = [1 — Tr a p 2 ] (p a 
is the reduced density operator of the a-th subsystem a = 0, ..N, where all the other degrees of freedom have been 
traced out) giving for the field (a = 0) and the fc-th ion entanglement the expressions below 

S L (po(t)) = 1 - ^P„ c (ao)^(ao)e- E "- 2| ^ |2(1 - cos ^ (m -" ) * ]) (26) 

n,77i 



S L ( Pk (t)) = 1 - ^^(a )P,;(ao) e ~ 2|Qfc|2(1 - cos[Afc(m -" ) * 1) , (27) 

where P°(a ) = e -l a °l 2 ^P and k = 1, ..N. 

We also consider the global entanglement measure in the form defined by Oliveira et al. [13] which for a system 
with M identical parties is given by 

EG ^=(dh)^J2 S M (28) 
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where Sl,(Po,) is the linear entropy of the a-th individual subsystem and d is the dimensionality of the each subspace. 
In our case, M = N + 1 and all subsystems are harmonic oscillators (d = oo) leading to 



= 1 ^ P n( a o) P ™( a o) | ^ g -2|a 3 | 2 (l-cos[A 3 (m-n)t]) _ p - £f =1 2|a, | 2 (l-cos[A 3 (m-n)*]) | _ ^ 



In Figfl] we plot the case with equal interaction strength of the ions with the field, otk = cto — lj an d = 10. In this 
equal coupling case, all the entanglements are very similar, with the field one slightly higher than the individual ion's 
(from Eq. (f2"T)) all ions have identical behavior). The global entanglement behaves like the individual entanglements. 
In Fig^a), we show that the field remains at a plateau of entanglement for longer and longer times as we increase 



0.8 




FIG. 1: Plots of entanglements as a function of time for all coherent initial state: global (continuous), field (dotted lines), and 
individual ion (dot-dashed lines) ; for the following parameter choices N = 10, Afc = l.OAi (fc = 2, .., 10) , on = oto = 1.0. 



the number of ions (N — 2, 10, 50), except at those times at which it disentangles or partially disentangles. On the 
other hand, the global entanglement changes little with N as shown in Figl^b). In Fig[3l we plot the case of unequal 
couplings with N = 5, where the incommensurability of the coupling constants causes the global entanglement curve 
to become irregular and non-periodic. It is interesting to notice that in this case the field remains in the plateau and 
does not return to zero. This happens because the field is always entangled with some of the ions, although individual 
ions present periodic behavior (entanglement of ion-2 is shown in this case) . 

B) Initially squeezed vacuum state for the Held: let us consider the initial state 1^(0)) = ®" =1 |aj) ® |0, £/), with 
£/ = Tf£ %vi ■ For this initial state we get a similar expression for the linear entropies of the field (|26[) and fc-th ion 
(f27|) . with the summation done only for even values of n, m and the distribution function P£ substituted by the one 
for squeezed state 



tanh(r^) n 
2 

(n even), where we have chosen <ff = 0. Therefore the GE in this case is also given by (f2l)|) with the same substitutions. 
In FigHl we plot the ion-2, field and global entanglements for N — 2 and unequal couplings with the field. The 
squeezing of the initial field state lowers the maximum values of the entanglement as compared with the coherent field 
state (see Figf2Ja)). Also, due to the fact that the couplings with the field are again different, the field and global 
entanglement exhibits non-periodic behavior as noticed in the previous cases (see FigEJ) ■ This allows one to engineer 
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(a) 

0.8 | . . r— 





A,,t 

FIG. 2: Entanglements as a function of time for all coherent initial state, (a) field, (b) global; for the following parameter 
choices N = 2, 10, 50, \ k = l.OAi (2 = 1, ..N) , an = a = 1.0. 

a decoherence scenario via cross-Kerr type interaction [33] of the field by choosing appropriate values of the parameter 
Xi = 2i] 2 gi/Ai for each ion. The easiest way of doing this modulation is by choosing the detunings appropriately 
through the application of external electric fields on each ion (static Stark effect) . 

C) Initially squeezed vacuum states for the ions: consider now all the ions initially squeezed |^(0)) = <8>™ = i|0,£j) <g> 
| a/) with £j = rje lVj (for simplicity ipj — for all j). The field and ion linear entropies result as 

^( P/ (t)) = i-nf =1 ^(^)^(^> 2|Qo|2(cos[A3(m " n)tH1) , (30) 

even n,m 

S L { Pk {t)) = l- ^(r fc )P^(r fc )e 2 l a °l 2 ( cos [ A ^ m -")*]-i). (31) 

even n,m 

The corresponding GE is given by the following expression 

1 N 
E G (1) = 1 - jj—j M^irMr^ + YP^rMr^e 2 ^ (-[^(m-n)*]-!) (32) 

even n,m j — 1 
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FIG. 3: Entanglements as a function of time for the all coherent initial state: global (dashed), field (dot-dashed lines) and 
ion-2 (continuous line); for the following parameter choices N = 5, A2 = \/2\i,\3 = l/v / 2Ai,A4 = \/3Ai,A5 = l/v3Ai, 
on = a = 1.0. 




FIG. 4: Plots of entanglements as a function of time for initially squeezed fieM and coherent ions: global (dashed), field 
(dot-dashed lines), ion-2 (continuous line),; for the following parameter choices N = 2, A2 = 1.25Ai , on — r c — 1.0. 



In Fig. ((5]) we show the time evolution of the linear entropy of the ion-2, field and the global entanglement. It is 
notable that the field and individual ion entanglements exchange the role from FigH] to FigEl namely that the field 
entanglement is now smaller than the entanglements of the individual ions. Also, it becomes clear that in the case 
where either the field or the ions are in the squeezed states, the subsystem which is coherent shows a typical periodic 
behavior, whereas the ones which are in the squeezed states have a more irregular behavior. 



V. CONCLUSION 



We have investigated the dynamics of a cluster of N trapped ions in the dispersive interaction with a quantized 
electromagnetic field. Depending on an appropriate preparation of the internal state of the ions in the cluster, we 
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FIG. 5: Plots of entanglements as a function of time for initially squeezed ions and coherent field: global (dashed), field 
(dot-dashed lines), and ion-2 (continuous line), for the following parameter choices N = 2, A2 = 1.25Ai , a = n = ri = 1.0. 



were able to obtain either cross-Kerr type interactions between field and ions or ionic dipole-dipole couplings that 
depend on the intensity of the motion of the trapped ions. With such dipole-dipole coupling, simple postselcction, 
and appropriate choice of the interaction time, one is able to generate non-local entangled mesoscopic states (NLMS) 
of two well separated ions, as we have shown. 

A simple preparation of all ions in the ground state has led to motion-dependent Stark-shifts in the form of cross-Kerr 
type interactions between field and ions. This allowed us to study both the bipartite entanglement of the field with 
the cluster or each of the ions with the rest of the system, and also to evaluate the multipartite global entanglement. 
We have considered coherent and squeezed initial states which are of interest in many physical situations and seems 
to be practically feasible in the current status of laser manipulated trapped ions experiments. In the case of a large 
cluster of ions, we analyze the possibility of observing controlled decoherence of the field in laboratory by choosing 
both an appropriate initial set of states for the system and an adequate set of effective interactions via detunings to 
dephase the periodicity of the entanglement of individual ions. 
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